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Abstract 

Using an algorithm for simulating equilibrium configurations, we study a fluctuating helical 
polymer either (i) contained in a cylindrical pore or (ii) wound around a cylindrical rod. We work 
in the regime where both the contour length and the persistence length of the helical polymer 
are much larger than the diameter of the cylinder. In case (i) we calculate the free energy of 
confinement and interpret it in terms of a worm-like chain in a pore with an effective diameter 
that depends on the parameters of the helix. In case (ii) we consider the possibility that one end 
of the helical polymer escapes from the rod and wanders away. The average numbers of turns at 
which the helix escapes or intersects the rod are measured in the simulations, as a function of the 
pitch pq. The behavior for large and small po is explained with simple scaling arguments. 

PACS numbers: 87.15.Aa, 36.20.Ey, 61.25.Hq, 82.70.-y 
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I. INTRODUCTION 



In this paper we study some of the equilibrium statistical properties of a confined helical 
or ribbon-like polymer. The cases of (i) a polymer contained in a cylindrical pore and (ii) 
a polymer wound around a cylindrical rod are considered. Some motivation is provided by 
the following observations: 

Biological polymers differ from synthetic polymers in that they are semi-flexible, with a 
persistence length much larger than the monomer size, and usually have a helical structure. 
This is well known for DNA, but F-actin also has a double-helical structure, while micro- 
tubuli are helical cylinders. The diameter of these biological polymers is in the range of 1 
to 25 nm. Polymeric helical structures are also found in self-assembling systems, consisting 
of either amphiphiles or peptides. In some cases the diameters and pitch lengths are much 
larger than for the biopolymers mentioned above. 

In amphiphilic systems the formation of helical ribbons has been observed in multicom- 
ponent mixtures of a bile salt or some other nonionic surfactant, phosphatidylcholine or 
a fatty acid, and a steroid analog of cholesterol The ribbons have typical diame- 

ters in the range of 5 to 20 /zm, and pitch angles between 10 and 60°. Other examples 
are ethanolic/ water solutions of diacetylenic phospholipids, in which the formation of hol- 
low tubules of diameter 0.6 /im and typical lengths of 10 to 100 /im has been observed 
H U, 0, 0] . Helically-coiled phospholipid-bilayer ribbons appear as metastable intermediates 
in the growth of these tubules. 

Other systems, which show spontaneous assembly of ribbons, are aqueous solutions of 
peptides (3, 0, 0, [l(3]. Depending on the solution conditions, the same peptide exists in 
different conformations, such as random coils, a-helices, or /3-sheets. At not too low peptide 
concentrations, the molecules self-assemble into long /3-sheet structures which form twisted 
ribbons (with a straight central axis). The width of these ribbons is about 4 nm, and their 
length is of the order of 500 nm P, Q. These ribbons can aggregate due to face-to-face 
attraction into twisted fibrils of a thickness of 8-10 nm. 

Interestingly, in a self-assembling system of gemini surfactants (two surfactant molecules 
covalently linked at their charged head group), the degree of twist and the pitch of the 
micrometer-scale ribbons has been found to be tunable by the introduction of opposite- 
handed chiral counterions 
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The confinement of polymers in cylindrical tubes is one of the classical problems in poly- 
mer physics. For biological polymers, such a confinement occurs, for example, when viral 
DNA of a bacteriophage squeezes through the narrow tail during DNA injection. Techno- 

H'cal advances in the manipulation of single molecules in micro- and nanofluidic devices 
I13I ] has fueled interest in the structure and dynamics of biological polymers in confined 
geometries 

Helical and twisted ribbons can be confined not only by external walls, but also by winding 
around each other, as in the fibril formation of twisted /3-sheet peptides mentioned above. 
The simple model we consider, consisting of a helical wound around a thin cylinder, is a step 
in this direction but leaves out some important physical features, such as the face-to-face 
attraction in the fibrils. 

II. FREE ENERGY OF CONFINEMENT 

The free energy AF of confinement of a fluctuating polymer of contour length £ in a 
cylindrical pore of diameter D is defined by 



Here Z(D,£) and Z(oo,£) are the partition functions of the polymer with one end fixed in 
the presence and absence of the cylindrical confining geometry, respectively. The quantity 
AF represents the work required to squeeze the polymer reversibly into the cylindrical pore. 
It may be evaluated in simulations by generating polymer configurations with one fixed end 
in an infinite volume with the Boltzmann probability, computing the fraction p(£) of the 
configurations of arc length I which lie entirely within a cylindrical domain of diameter D, 
and making use of Eq. (0) 

For a flexible, self-avoiding polymer with vanishing bending rigidity, AF is purely en- 
tropic. The confinement of such a polymer in a cylindrical pore is considered in Refs. 



In the worm-like chain model of a semi-flexible polymer, the polymer is represented by 
an inextensible line or filament r(s) with contour length £ and elastic energy 



exp(-AF/k B T) 



Z(D,£) 
Z(oo,£) 



p(£) 



(1) 





(2) 
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Here s specifies distance along the contour, t3 = dr/ds is the unit tangent vector, k is the 
bending rigidity, and P = n/ksT is the the persistence length. In the narrow-pore, long- 
polymer limit D ^ P ^ £, the polymer is almost a straight line, i.e. the angle between the 
tangent vector t 3 and the z axis or symmetry axis of the cylinder is a small quantity. In 
this case the right-hand side of Eq. (JTJ) decays as 

p(£) ~ e~ E ° l (3) 

for large £, where exp(— E^dz) is the largest eigenvalue of the transfer matrix of a slice of 
the system of thickness dz. The quantity Eq 1 represents a typical contour length at which 
the configurations intersect the pore wall. According to Eqs. (0) and (jBJ) the confinement 
free energy per unit length Af = AF/£ is given by 

A/ E (P,D) = -^, (4) 



k B T uv ' ' pwdw 

where the dependence on P and D follows from simple scaling or dimensional arguments 
jisl llil 2f| . Similarly, for a pore with a rectangular cross section with edges L\, L 2 <C P, 

A/ Wi,i 2 ) = 4|-^ -^r! • <- r >) 




k B T pl/3^2/3 

The quantities A Q and An on the right-hand sides of Eqs. (J1J and © are dimensionless 
universal numbers A Q and A n , which are the same for all worm-like chains. 

The prediction A n = 1.1036 was obtained in Ref. jscj by solving an integral equation 
numerically which arises in an exact analytical approach. Measuring the probability p(£) 
in Eq. Q in simulations, fitting the large £ behavior with the exponential form Q, and 
making use of Eqs. and Bicout and Burkhardt 2l] estimated 

A = 2.375 ±0.013 , A n = 1.108 ±0.013 . (6) 

An earlier estimate from simulations, A a = 2.46 ± 0.07, was given by Dijkstra et al. Q|. 



III. HELICAL POLYMER MODEL 



In this paper we generalize the above results to helical polymers or chiral ribbons, which 
have spontaneous curvature and torsion. Again the polymer is replaced by a curve r(s) 
of fixed contour length S. To each point on the line a right-handed triad of unit vectors 
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ti(s), t 2 (s), t 3 (s) is assigned, where t 3 = dr/ds is the tangent vector and ti, t 2 correspond 
to principal axes of the polymer cross section. The rotation of the triad along the curve is 



governed by the generalized Frenet equations 

dti 



22, 



24J 



ds 



WXtj, OJ — txUi + t 2 UJ 2 + t 3 u 3 



or 



The elastic energy is given by J221, 13, |3] 



(7) 



(8) 



E 



helix 



1 y 3 r r S 

j=l J o 



(9) 



where the coefficient b\ and b 2 are bending rigidities along the principal axes of the cross 
section, and 63 is the twist rigidity. The parameters Uj(s) and uj j(s) determine the curva- 
tures and torsions in the deformed and stress-free states of the polymer, respectively. Since 
the energy is quadratic in the deviations 5ujj = ujj — the distribution of Sujj is Gaussian, 
with zero mean and second moment 



(5ui(s)6uj(s')) 



Sij S(s - s) . 



(10) 



We restrict our attention to the case u>oj(s) = constant, corresponding to a helical polymer 
with spontaneous curvature and torsion but without spontaneous twist. In the absence of 
fluctuations, i.e. in the limit b\ = b 2 = b 3 = 00, the Frenet equations are readily solved 
yielding 



r(s) = r(0) + — { t 3 (0)sin(cc>os) + e(0) w 03 



+ e(0) x t 3 (0) [1 - cosmos)] } , 



sin (u s) 



(11) 



where 



e(s) = ti(s) h t 2 (s) h t 3 (s) , uj = Hi + ^02 + ^03) 

LOO LOq CUq 



, w 02 



, <^03 



2 \l/2 



;i2) 



Equation (fTTj) represents a helix with radius r and pitch p , where 



winding around an axis pointing in the direction of the unit vector e(0). 

Including Gaussian fluctuations according to Eq. (JTUJ), Panyukov and Rabin 0, Q 
showed that 

(t 4 ( S )-t J (0)) = ( e - rs ) y . , (14) 
where T is the matrix with elements 



V / k 



k^ok ■ (15) 

The two-point correlation function of the unit vector e(s) in Eq. ()12|). which is directed 

along the axis of the helix, follows from this result. In the special case b = b\ = 62 = &3 
considered in our simulations, 

(e( S ) • e(0)> = e'^ , L p = ^ , (16) 

where L p is the persistence length. 



IV. SIMULATIONS 



24j, we replace the differential equations (jHJ) by the difference 



Following Kats et al. 
equations 

t ik (s + ds) = Oijt jk (s) (17) 
3 

in our simulations. Here t ik denotes the fc-th component of tj with respect to a fixed Cartesian 
coordinate system, O is the orthogonal matrix 

2 J V 2 



0= [1 + -Ads] (l--Ads) , (18) 



and A is the antisymmetric matrix with elements Aij = tijkUJk- The difference equations 
are consistent with the Frenet equations (JHJ to first order in ds, and the orthogonality of 
the matrix O preserves the orthonormality of the tj in the simulations. 

For simplicity we set b = b\ = 62 = bs, corresponding to Eq. f)16|) . In accordance with 
Eq. (JTUJ), the 5uij(s) are chosen randomly from a Gaussian distribution with zero mean and 
standard deviation (fceT ' /bds) 1 / 2 , where ds <C L p = b/ksT. 
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V. HELICAL POLYMER IN A CYLINDRICAL PORE 



We have determined the confinement free energy of a helical polymer fluctuating in a nar- 
row cylindrical pore from simulations. Cylinders with both circular and square cross sections 
were considered, and we use the same symbol D for the diameter and edge, respectively 
The symmetry axis of the cylinder defines the z axis of our fixed coordinate system. 

The helical polymer was generated step by step using the numerical procedure described 
in the preceding Section. The radius r and pitch p were set to desired values by choosing 
u i, u 03 in Eq. (JT5j) appropriately, with u 02 = 0. The starting point was chosen randomly, 
apart from the requirements that the stress-free helix fit inside the cylinder, with its axis 
parallel to the z axis. This is the case for the initial vectors ti(0) = (uj 03 /uj ,0,lj i/uj ), 
t2(0) = (0,1,0), and t 3 (0) = (—^01/^0,0,^03/^0), which were used. The other simulation 
parameters were ds = 10~ 4 , D = 1, L p — bi/k B T = b 2 /k B T = b 3 /k B T = 8000. Clearly 
ds <C D <C L p . The pore is narrow in comparison with the persistence length, and ds is 
small in order to approximate the continuum model (|7jl. 

To obtain the free energy of confinement of the helical polymer, we proceeded as discussed 
below Eq. (JTJ, generating many polymer configurations and computing the probability P(n) 
that the polymer has not yet intersected the pore wall 26| after n steps of the algorithm. The 
determination of P(n) was based on 50, 000 independent helices. For large n an exponential 
decay 

P{n) ~ e~ Xon , (19) 

similar to the result Q for semi-flexible polymers, is expected. According to Eq. ([TJ, the 
free energy of confinement per unit length along the axis of the helix is given by 

*L = A £ = Po f2 0) 

k B T ids' ? [pl+(2nr y] 1/2 ' 1 ' 

Here we use the relation i = £s between the contour length s and the corresponding length 
I along the axis of the helix. The persistence length L p , defined with respect to the contour 
length as in Eq. pfijl . and the corresponding persistence length P, defined with respect to 
length along the axis of the helix, also satisfy P = £L P . 

A helical polymer with persistence length L p in a pore with diameter D -C L p has the 
same confinement free energy as a semi-flexible polymer with persistence length P = £L P in 
a pore with effective diameter D eS . To define D cfi quantitatively, we equate the free energies 



of confinement (J3J) and (j20j). obtaining 



4, _(^ P ) 1/3 X A D (^ P )V\ 



£ 2 / 3 fds u ' D 2 / 3 2£ds 

The probability P(n) for ro = po = 0.3, with the other simulation parameters noted 
above, is shown in Fig. 1. The data are in good agreement with the exponential decay ()19j) . 
and the values of Aq are given in the figure caption. As in the case of a semi-flexible polymer 



2l| . the curves P(n) for the circular and square cylinders practically coincide when plotted 
versus Xon instead of n. 

For the exponential decay ()19|). the mean number of steps of the algorithm at which the 
polymer intersects the wall equals Ag , corresponding to Ni = £ds/(po\o) turns of the helix. 
The values of Ao in the caption of Fig. 1 yield iV» = 8.0 and 8.7 for the circular and square 
cross sections. Since the number of turns before intersecting the wall is fairly large, the helix 
should be equivalent to a worm-like chain in a pore of width D c d = D — 2r . To check the 
equivalence quantitatively, we use Eq. (J2TJ) with D e g = D — 2r and the values of Ao in the 
caption of Fig. 1 to predict the amplitudes A Q , A n . This yields 

A = 2.45 ±0.05 , A n = 1.12 ±0.04 , (22) 

in good agreement with the results (JBJ) for semi-flexible polymers. 

We have also studied the dependence of D cS on the polymer pitch p , keeping the radius 
ro and the persistence length L p constant. For small p$ the polymer makes many turns before 
intersecting the wall and is equivalent to a semi-flexible polymer in a pore of diameter D— 2ro, 
as discussed in the preceding paragraph. In the limit p Q — > oo, the helical polymer does not 
make any turns before intersecting the wall and corresponds to a semi-flexible polymer in a 
pore of diameter D. As po increases from to oo, D c r is expected to vary monotonically 
between these two limiting values. 

For various values of po we have computed the probability P(n) that a helical polymer 
with radius ro = 0.3, persistence length L p = 8000, and contour length nds in a cylindrical 
pore with a circular cross section of diameter D = 1 does not intersect the wall. The 
corresponding Ao was obtained from an exponential fit for large n. Finally _D e fr was 
calculated using Eq. (|2~Tj) and the best estimate (JBJ) for A . The results are shown in Fig. 2. 
The data do indeed interpolate between the expected limiting values D — 2r = 0.4 and 
D = 1 for small and large p , respectively. 
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The crossover region in Fig. 2, where D c g varies most rapidly with p$, is around po ~ 40, 
D e g « 0.7. According to Eqs. (j3J), (jlj), and pPjl. these values of po and D c fr correspond to 
Ni = (^Lp) 1 / 3 -Dpff 3 / (AoPo) ~ 0.2 turns of the helix before intersecting the wall. 

VI. HELICAL POLYMER ENCIRCLING A CYLINDRICAL ROD 

In this Section we consider a helical polymer wound around a long cylindrical rod with 
a circular cross section and diameter D ^ L p . We study the possibility that the polymer 
generated in the simulation escapes from the rod as n increases and wanders away. 

In the simulations the parameters ds = 10~ 4 , bi/ksT = b^lksT = b^/ksT = L p = 8000 
were the same as in the preceding section. The diameter of the rod was D = 0.2, and the 
radius of the helix was r = 0.3. For these parameters ds D < 2r <C L p . The starting 
point of the polymer was chosen randomly, apart from the requirements that the stress-free 
helix wind around the cylindrical rod without touching it, with the axis of the helix parallel 
to the rod. 

^From the simulation data we computed the probability P(n) that after n steps the 
polymer has not yet intersected the rod |26[. Each curve P(n) is based on 10, 000 independent 
helices. The results for three different values of the pitch p are shown in Fig. 3. Unlike 
the case of a polymer in a cylindrical pore, shown in Fig. 1, P(n) does not decay to zero 
as n increases. Instead, above a characteristic value which depends on the pitch, the curve 
flattens and approaches a nonzero limiting value. This is because the polymer generated in 
the simulation sometimes escapes from the rod, due to a sufficiently large fluctuation, and 
wanders away as n increases, without ever returning to intersect the rod. 

A simple theory of the escape, which suggests P(n) = A + Be~ Cn , in qualitative consis- 
tency with Fig. 3, is given in the Appendix. 

We determined the average number of turns at which the helix escapes from the rod or 
intersects it by making two checks after each step of the growth algorithm: (i) If the distance 
of the endpoint r(s) from the axis of the rod is less than D/2, the polymer has intersected 
the rod. (ii) If the distance of the endpoint r axis (s) of the axis of the helix is greater than 
ro + -D/2, the circular cross section of the helix no longer encircles the rod, i.e. the helix 
has escaped. Geometrically r axis (s) is determined as follows: Since the unit vectors ts(s) 
and e(s) are tangent to the helix and directed along its axis, respectively, e(s) x ts(s) is 
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directed perpendicularly from the point r(s) on the helix contour toward the corresponding 
point r axis (s) on the axis of the helix. Thus, 

e(s) x t 3 (s) 



r(s) +r 



e(s) x U(s) 



= r (a)+ti(a)^-t 2 ( a )^, (23) 

where we have used Eqs. (fT2j) and (fT3*j). 

In Fig. 4 the average numbers of turns AT*, iV| at which the helix escapes from the 
rod and the average number of turns Ni at which the helix intersects the rod are shown 
as functions of p . For each value of p , 10,000 independent configurations were generated. 
Each configuration was continued until it intersected the rod |26J] or the number of steps 
of the algorithm exceeded 5 x 10 6 , whichever came first. The average iV* is based on all 
configurations which escape, independent of whether they return to intersect the rod or not. 
The quantity N 2 is the average value for only those configurations which escape and in 
5 x 10 6 steps of the algorithm still have not intersected the rod. 

For p Q < 1 the probability of the polymer escaping from the rod is so small that N]' 2 
could not be determined reliably with configurations of 5 x 10 6 steps. For larger p , the data 
for iVg 1 and N 2 practically coincide, indicating that the polymer rarely returns to intersect 
the cylinder once it has escaped. The data are in excellent agreement with N^' 2 , iVj ~ p$ 1 
for large po and Ni ~ p for small pq. These power laws may be understood as follows. 

The transverse fluctuations of the endpoint of the axis of an unconfined helical polymer 
of length £ and persistence length P about the corresponding endpoint of the unstressed 
helix are readily calculated from Eq. ()16|) and given by 27] 

(rl) = \ l -p (24) 

for I < ?. Here both I = £s and P = £L p are measured along the axis of the helix, as 
discussed below Eq. (|20p. Equation ()24j) also applies to the worm-like chain. Apart from the 
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IS, 2fl. The 



factor 2/3, Eq. ()24|) follows from simple scaling or dimensional arguments 
powers of r_|_, £, and P in Eq. (J24|) are also consistent with E £ ~ 1 and r± ~ D in Eqs. (J3J) 
and flU). 

Replacing r± in Eq. (|24|) by ro + D/2 = 0.4, as in the simulations, and solving for 
N e = £/p yields the estimate 

12 

p ' {pi + 3.6) 
10 



for the number of turns of the helix at which the typical transverse displacement equals the 
value needed for escape from the rod. Roughly speaking, the polymer wrapped around the 
rod escapes in N e turns, as given by Eq. (J2o"|) for N e ~ 1, i.e. p ~ 12. For smaller p , the 
typical transverse fluctuations in a single turn of the helix are too small for escape from the 
rod, and the helix is more likely to intersect the rod than to escape. Equation (|25|). which 
ignores this possibility, no longer applies. As noted above, for p < 1 the escape probability 
is too small for a reliable determination of N 1 ' 2 with configurations of 5 x 10 6 steps. 

In the region N e ~ 1, i.e. Pq ~ 12, Eq. (J23j) . which corresponds to the solid curve in 
Fig. 4, is in good quantitative agreement with the simulation data for A e 1,2 . For large po, 
N e pa 12/p Q . The coefficient 12 is an order-of-magnitude estimate that happens to give a 
good fit to the simulation data, whereas the power law N e ~ p^ 1 for large po is exact. An 
argument based on Eq. (j2*4^) similar to the one for N e predicts N ~ p$ 1 for large po, in 
excellent agreement with Fig. 4. Note that the data points for A e 1,2 and iVj practically 
coincide. 

Since the possibility of escape is negligible for small po, the helical polymer is equivalent 
to a worm-like chain in a pore with diameter D e g = 2tq — D. To estimate the average number 
of turns Ni at which the polymer intersects the rod, we replace D by _D e g in Eqs. Q and 
(J3J and solve for the typical intersection length £ pa Eq 1 . This yields i ~ P 1 / 3 . Substituting 
r± ~ D c q in Eq. an d solving for i leads to the same result. According to the discussion 
below Eq. (J2DJ), P = ^L p pa L p (p /2iir ) for p r . Keeping track of the powers of po, 
we obtain the power law Ni = £/po ~ Pq 2 ^ f° r small po, in excellent agreement with the 
simulation data in Fig. 4 

In Fig. 5 the fractions f 2 , and /, of the 10,000 configurations which contribute to 
iVg 1 , N 2 , and Ni in Fig. 4 are shown as functions of Pq. For po ~ 1, / e 1,2 pa and fi ~ 1, i.e. 
almost all the configurations intersect the cylinder and never escape. Around po pa 10 the 
curves cross, and for larger po the polymer is more likely to escape than to intersect the rod. 

VII. CONCLUDING REMARKS 

We have studied some statistical properties of a helical polymer in cylindrical restrictive 
geometries of diameter D, in the limit that the persistence length P along the axis of the 
helix is large in comparison with D and the radius tq of the helix. In this limit the helical 
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polymer has much in common with the worm-like chain. We interpret the simulation data 
for the free energy of confinement in a cylindrical pore using the scaling form (£Q) for a 
worm-like chain in a pore, with an effective diameter D e g that is renormalized by the helical 
structure. As the pitch p of the helix increases from to oo, -D e fr increases monotonically 
from D — 2r to D, as shown in Fig. 2. 

Thinking in terms of a worm-like chain also proves useful in connection with the escape 
of the helical polymer encircling a cylindrical rod. In the limit P ^> tq the transverse 
fluctuations of the axis of the helix are given by the same result (|24|) as for the worm-like 
chain. As po increases, the typical transverse displacement in one turn of the helix also 
increases, resulting in a greater probability per turn of escape. We have used Eq. (J24*Jl to 
estimate the average number of turns at which the helix escapes from the rod or intersects 
it. The simulation data in Fig. 4 are in excellent agreement with the predicted asymptotic 
forms N e , Ni ~ p$ 1 , Aj ~ p 2 ^ 3 for large and small p , respectively. 

It would be interesting to include an attractive interaction between the rod and the poly- 
mer wound around it. In the fibril formation mentioned in the Introduction, the attraction 
is an essential ingredient. 
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APPENDIX A: SIMPLE THEORY OF ESCAPE OF A POLYMER 

Let us define P^ e as the probability that the polymer has neither intersected the cylindri- 
cal surface of the rod nor escaped in the first N turns of the helix and Pfj as the probability 
that it has not yet intersected the cylindrical surface but that that it has escaped. Treating 
each turn of the helix as statistically independent, we denote the probability that a polymer 
which has not yet escaped does escape in the next turn by q, and the probability that it 
neither escapes nor intersects the rod in the next turn by p. (The third possibility, that it 
intersects the rod in the next turn, has probability 1-q-p.) In addition we assume that once 
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the polymer escapes, it never intersects the rod. 
These assumptions imply the recurrence relations 



PnXi = pPn , (Al) 
P e N+l = qPST + Pn , (A2) 

with initial conditions P^ e = 1, Pq = 0. Writing down the first few iterates, it is easy to see 
that 

P N e = P N , (A3) 
Pn = <1 T Z^- (A4) 

The probability P^ = P^f + Pfj that the polymer has not yet intersected the rod after iV 
steps is analogous to P{n) in Section V. From Eqs. (|A3jl and (|A4jl 

Pjv = ^ + — ; P • (A5) 

1 — p 1 — p 

Thus, as N increases, Pn decays exponentially from Pq — 1 to Pqo = q/(l —p)- The mean 
number of turns N e at which escape occurs is given by 

An analogous calculation for the mean number of turns Ni at which the polymer insects the 
rod yields Ni = N e . 

This theory is obviously an oversimplification, but the form (|A5|) of the decay, P/v = 
A + Be~ CN , is qualitatively consistent with Fig. 3. 
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FIG. 1: Probability P(n) that a helical polymer with radius ro = 0.3, pitch po = 0.3, and persistence 
length L p = 8000 in a cylindrical pore does not intersect the pore wall in the first n steps of the 
algorithm. The full circles (•) correspond to a pore with a circular cross section with diameter 
D = 1 and the triangles (A) to a square cross section with edge length D = 1. Fitting the data to 
Eq. (|19j) for large n yields Ao = 6.57 x 10 -6 and 6.01 x 10 -6 , respectively. The full line corresponds 
to the exact exponential decay e _A ° n . 
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FIG. 2: The effective diameter D e s as a function of the pitch po for a helical polymer with radius 
ro = 0.3 and persistence length L p = 8000 in a cylindrical pore with a circular cross section with 
diameter D = 1. The data interpolate between the limiting values D — 2r$ = 0.4 and D = 1 for 
small and large po, respectively. 
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FIG. 3: Probability P(n) that a helical polymer with radius = 0.3 and persistence length 
L p = 8000, wound at the fixed end around a cylindrical rod of diameter D = 0.2 does not intersect 
the rod in the first n steps of the algorithm. The pitch of the helix is po = 10 (•), 30 (A), and 100 

(°)- 
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FIG. 4: Average numbers of turns (A), (•) at which the helix escapes from the rod, and the 
average number of turns iVj (o) at which the helix intersects the rod as a function of po- Here is 
based on all the configurations which escape, independent of whether they return to intersect the 
rod or not; N% is based on the configurations which escape and in 5 x 10 6 steps of the algorithm 
do not return to intersect the rod. The dashed lines on the left and right have slopes -2/3 and -1, 
respectively, and the solid line shows the prediction 1)25)1. 
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FIG. 5: Fractions f} (A), /| (•), and fi (o) of the 10,000 configurations which contribute to AT*, 
N% , and AT, in Fig. 4, as a function of p®. 
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